Electron spin relaxation of NQCeo in CS2 



o 
o 

(N 

a 

in 



C/2 



John J. L. Morton,^' ^ 'El Alexei M. Tyryshkin,^ Arzhang Ardavan,^ 

Kyriakos Porfyrakis,^ S. A. Lyon,'^ and G. Andrew D. Briggs^ 

'Department of Materials, Oxford University, Oxford 0X1 3PH, United Kingdom 
'^Clarendon Laboratory, Department of Physics, Oxford University, Oxford 0X1 3PU, United Kingdom 
''Department of Electrical Engineering, Princeton University, Princeton, NJ 08544, USA 

(Dated: February 2, 2008) 

We examine the temperature dependence of the electron spin relaxation times of the molecules 
N@C6o and NQCto (which comprise atomic nitrogen trapped within a carbon cage) in liquid CS2 so- 
lution. The results are inconsistent with the fluctuating zero field splitting (ZFS) mechanism, which 
is commonly invoked to explain electron spin relaxation for S > 1 spins in liquid solution, and is the 
mechanism postulated in the literature for these systems. Instead, we find an Arrhenius temperature 
dependence for NQCeo, indicating the spin relaxation is driven primarily by an Orbach process. For 
the asymmetric NQCto molecule, which has a permanent ZFS, we resolve an additional relaxation 
mechanism caused by the rapid reorientation of its ZFS. We also report the longest coherence time 
(T2) ever observed for a molecular electron spin, being 0.25 ms at 170K. 

PACS numbers: 76.30.-v,81.05.Tp 
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INTRODUCTION 

The encapsulation of atomic nitrogen within a fuUerene 
shield has provided a uniquely robust molecular electron 
spin 0. Its unique relaxation properties have enabled 
the observations of a novel type of electron spin echo en- 
velope modulation (ESEEM) and attracted attention 
as a potential embodiment of a bit of quantum informa- 
tion Q. 

In high spin systems {S > 1) in liquid solution, a 
fluctuating zero field splitting (ZFS) has habitually been 
cited as the dominant relaxation mechanism since tran- 
sition metal ions were first studied by EPR 0, 0| . When 
relaxation in NQCgo (which has electron spin S — 3/2) 
was first studied, it was therefore natural to assume that 
the same ZFS mechanism applied 0. However, to date 
there has been little evidence to support this hypothe- 
sis. For example, no temperature dependence has been 
reported for NCSCeo in solution; such a study is critical 
in determining unambiguously which relaxation mecha- 
nisms are relevant. Measurements have been reported 
in CS2 and toluene solutions 0; however, the analysis 
of these results ignored the effects of magnetic nuclei in 
toluene, which we have found to contribute significantly 
to the relaxation . Finally, the previous measurements 
were performed using fuUcrene solutions that were suffi- 
ciently concentrated for (C6o)n aggregates to form, so it 
is difficult to conclude which phase (liquid or solid) the re- 
ported T1/T2 times correspond to Consequently, the 
favoured relaxation model of a zero-field splitting (ZFS) 
fluctuation has little direct evidence to support it, and 
must be critically re-evaluated. 

In this letter we report relaxation times for both 
N@C6o and NQCyo in CS2 solution, which, conveniently, 
lacks nuclear spins in the dominant isotopes of its con- 
stituents. We find that the temperature dependence of 
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FIG. 1: Continuous wave EPR spectrum of NQCeo in CS2 
at room temperature. Each line in the triplet signal is la- 
beled with the corresponding projection Mj of the nu- 
clear spin. Measurement parameters: microwave frequency, 
9.67 GHz; microwave power, 0.5 fiW; modulation amplitude, 
2 mG; modulation frequency, 1.6 kHz. 



the relaxation times is inconsistent with the previously 
proposed ZFS mechanism, and suggest an alternate Or- 
bach relaxation mechanism. We extract an energy gap 
which matches well the first excited vibrational state of 
the fuUerene cage. 



MATERIALS AND METHODS 

High-purity endohedral NOCgo was prepared 0, dis- 
solved in CS2 to a final fuUerene concentration of 1- 
2-10^^/cm'^, freeze-pumped in three cycles to remove oxy- 
gen, and finally sealed in a quartz EPR tube. The 
fuUerene concentration used (~ 1 fjM) was well below the 
cluster formation threshold [9j. Samples were 0.7-1.4 cm 
long, and contained approximately 5 • 10"'^^ NQCeo spins. 
Pulsed EPR measurements were performed using an X- 
band Bruker Elexsys580e spectrometer, equipped with 
a nitrogen-flow cryostat. T2 and Ti times were mea- 
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sured using 2-pulse (Hahn) electron spin echo (ESE) and 
inversion recovery experiments, respectively. The 7r/2 
and TT pulse durations were 56 and 112 ns respectively. 
Phase cycling was used to eliminate the contribution of 
unwanted free induction decay (FID) signals. 

Figured shows the continuous- wave EPR spectrum of 
NQCgo in CS2 at room temperature. The spectrum is 
centered on the electron g-factor g — 2.0036 and com- 
prises three narrow lines (linewidth < 0.3 /iT) resulting 
from the hyperfine coupling to "'^'^N The relevant 

isotropic spin Hamiltonian (in angular frequency units) 
is 

Ho ^ ujeSz ~ ujjlz + a-S-I, (1) 

where We = g(3B()/h and loj = giP„BQ/h are the electron 
and "'^'^N nuclear Zeeman frequencies, g and gi are the 
electron and nuclear g-factors, /3 and /3„ are the Bohr 
and nuclear magnetons, h is Planck's constant and Bq is 
the magnetic field applied along z-axis in the laboratory 
frame. Each hyperfine line (marked in Fig.^with Mj = 
and ±1) involves the three allowed electron spin transi- 
tions AMs = 1 within the S = 3/2 multiplct. These 
electron spin transitions remain degenerate for Mj = 
but split into three lines for Mj = ±1. This additional 
splitting of 0.9 fiT originates from the second order hy- 
perfine corrections and leads to a modulation of the elec- 
tron spin echo decay 

RELAXATION OF NOCeo IN CS2 

Spin relaxation times Ti and T2 for NQCgo in CS2, 
measured on the central Mj = hyperfine line, are shown 
on a logarithmic scale in Figure 13 for a range of tem- 
peratures (160K to 300K), demonstrating an exponen- 
tial temperature dependence and a roughly constant ratio 
T2 sa (2/3)Ti over the full temperature range. This con- 
trasts with previous findings which reported no tempera- 
ture dependence for T2 Below 160K, the CS2 solvent 
freezes as a polycrystal, leaving regions of high fuUerene 
concentration around grain boundaries. This dramati- 
cally increases the local spin concentration, and T2 be- 
comes extremely short due to dipolar s pin coupling (the 
so-called instantaneous diffusion effect [iM Ea. Il4l | ) . 

As this is an 5* = 3/2 spin system, one might expect 
several different relaxation times corresponding to the 
different AMs = 1 transitions. However, in the exper- 
iments presented in Figure |21 all decays were well de- 
scribed by monoexponentials. Given two similar expo- 
nential decays, it is notoriously difficult to extract any- 
thing other than a single, average decay constant from 
an exponential fit. Here, we take advantage of a recently 
reported mechanism for electron spin echo envelope mod- 
ulation (ESEEM) to separate the relaxation times for 
different electron transitions. This modulation generates 
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FIG. 2: Electron spin relaxation times (Ti and T2) of 
N@C6o in CS2, measured using the central Mj = line. The 
ratio T2 « (2/3)Ti is maintained over the full temperature 
range for which the solvent remains liquid. 

an echo intensity for transitions on the Mj = ±1 lines 
which varies as a function of the delay time, r, as 

t^M,=±i(T) = 2 + 3cos2(5t. (2) 

The oscillating component arises from the 'outer' coher- 
ences (from the Ms — ±3/2 : ±1/2 transitions), whilst 
the unmodulated component arises from the 'inner' co- 
herences (from the Ms = +1/2 : —1/2 transition). If 
T2 relaxation is included, Eq. |21 transforms to: 

VMr=±i{r) = 2exp(-2r/T2,,)+3exp(-2T/T2,o)cos2^r, 

(3) 

where T2,i and T2,o are the relaxation times of the 'inner' 
and 'outer' coherences, respectively. Thus, by fitting to 
the modulated ESEEM decay, the individual relaxation 
times T2,i and T2,o can be extracted. Ti and T2 times 
measured for the high- field (M/ = — 1) hyperfine line are 
shown in Figure|31 Ti was measured in the standard way 
(inversion recovery) , and so only one (average) value was 
obtained. 

The behaviour of Ti appears identical for both central 
and high-field lines, indicating that relaxation caused by 
the hyperfine interaction with the nitrogen nuclear spin is 
negligible. The T2,i measured on the high-field Mj = — 1 
hyperfine line correlates closely with the T2 measured 
on the central Mj = line. Remarkably, both of these 
T2 times remain approximately 2/3 of Ti over the full 
temperature range studied. For the high-field line, the ra- 
tio of T2,o to T2,i also stays constant at about 2/3. The 
fact that certain ratios between Ti, T2,i and T2,o remain 
constant over a broad temperature range is a strong in- 
dication that all of these relaxation times are limited by 
the same mechanism. In the following section, we review 
different relaxation mechanisms which might account for 
the observed temperature dependence. 
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FIG. 3: Electron spin relaxation times (Ti and T2) of 
N@C60 inCS2, measured using the high field AIi = —1 line. 
ESEEM is used to resolve the individual decay rates of the 
inner and outer coherences (see Eq. |HJ . Dashed curves show 
corresponding data taken for the central Mj = line, for 
comparison. 



ZFS fluctuations 



result from rotational tumbling in molecules that have a 
permanent non-zero ZFS (such as in N@C7o). In the case 
of a degenerate S = 3/2 system, a fluctuating ZFS term 
leads, in general, to two different decoherence times Ts^, 
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for the transitions that we refer to here as 'inner' and 
'outer' respectively. D^^^ = _|_ 3^2^ g^^^ ^ g^^.^ 
the coupling and rhombicity ZFS parameters, Tc is the 
correlation time of the fluctuations, and LUe is the electron 
spin transition frequency. 

The predicted Ti times arising from the same mecha- 
nism are: 



(Tm)" 



-D2 
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(6) 



Spin relaxation is manifested in fluctuating terms in 
the spin Hamiltonian and arises from fluctuating mag- 
netic dipoles (either nuclear or electronic), and other mo- 
tions causing variations in the interactions between the 
spin and its environment. The trapping of endohedral ni- 
trogen in a high symmetry environment suppresses most 
of the conventional spin relaxation mechanisms (zero- 
field splitting (ZFS) interaction, anisotropic g matrix, 
electron-nuclear dipolar coupling and nuclear quadrupole 
interaction). Indeed, it has been proposed that the 
dominant relaxation process arises from small deviations 
from this ideal symmetric environment, caused by cage 
deformations from collisions with solvent molecules 
For example, the modulation of the hyperfine interac- 
tion through such collisions is a possible relaxation path- 
way. This was dismissed in earlier reports on the basis 
that the expected Mj dependence of linewidth that this 
mechanism predicts is not observed ij. However, as all 
linewidths are likely to be instrument ally limited, this 
observation did not constitute a rigorous confutation. 

The mechanism favoured in the literature is that of 
a ZFS fluctuation, again caused by deformation of the 
spherical Ceo cage through solvent collisions 0. Given 
the concentrations of fuUerene solution that were re- 
ported in these earlier studies, a large amount of fullerene 
aggregation is expected 0] and so it is unlikely that the 
NQCgo molecules being studied had any direct contact 
with solvents. Nevertheless, deformations of the cage, 
through whichever mechanism (such as collisions with 
other Cgo molecules in the cluster) , will give rise to some 
time- varying ZFS. Alternatively, ZFS fluctuations may 
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The individual values of Ti^i and Ti.o cannot be re- 
solved in a simple inversion recovery experiment, and 
thus only their average can be determined (with respec- 
tive weights 2 and 3). 

In the fast tumbling limit {uJeTc << 1), the theory pre- 
dicts these two Ti times to be identical, and equal to both 
types of T2, contrary to our observed ratio of 2/3. Mov- 
ing away from the fast-tumbling limit, values for -De// 
and Tc can be derived given any values for Ti and T2. 
Since the ratio between these times is dictated purely by 
Tc, the fact that the ratios stay fixed implies Tc, the cor- 
relation time of the ZFS fiuctuations, stays fixed over the 
broad temperature range (160 to 300K). This would be 
surprising, as the viscosity of CS2 changes by an order 
of magnitude over this temperature range |l(il |. Thus, we 
conclude that the previously suggested ZFS fluctuation 
mechanism cannot explain the observed temperature de- 
pendence of Ti and T2, nor their mutual correlation, and 
therefore seek alternative explanations for the behaviour 
observed. 



Orbach relaxation process 

The temperature dependence of Ti is well described by 
an Orbach relaxation mechanism (see Figure 0}. This is 
a two-phonon relaxation process whose energies are reso- 
nant with a transition to an excited electronic state (i.e. a 
vibrational or orbital state which lies outside of the space 
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log(1 /(Ti +A)) + log(A) = - 375 / T 
A = 51 .5 
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FIG. 4: The temperature dependence of Ti of NQCeo is 
linear in Arrhenius coordinates, consistent with the Orbach 
relaxation mechanism. An energy gap A — 32(1) meV 
= 375K can be extracted. Because we cannot make a low- 
temperature approximation in this case, the standard Orbach 
plot of log(l/Ti) vs. 1/T must be adjusted to include the 
constant of proportionality, A (see Eq. |HJ . The plot is then 
recursively fit to fine-tune A and obtain the slope, A/k. Ti is 
given in microseconds. 



considered by the spin Hamiltonian) . The Ti tempera- 
ture dependence is dictated by the distribution of phonon 
energies, and is of the form: 



Ti 



.A/kT 



1) 



(8) 



where A is the energy gap to the excited state and A 
is some constant which involves terms associated with 
spin-orbit coupling (and therefore with the ZFS, ^"^N hy- 
perfine coupling and g-tensor in the excited state) 17] . A 
fit to the data in Figure H yields A = 32(1) meV. This is 
a close match to the energy of the first vibrational mode 
of Cgo (273 cm~^, or 34 meV) which has been theoret- 
ically calculated and observed b y R aman spectroscopy 
of Ceo in CS2 solution at 300K [II Ei^l, indicating 
that this may be a vibrational spin-orbit Orbach pro- 
cess 0,1131 • This first excited vibrational mode, termed 
Hg{l), breaks the spherical symmetry of the molecule, 
reducing it to axial. The small difference between A 
observed here compared with that seen in the Raman 
spectroscopy of Cgo could be due to a shift in vibrational 
energies due to the presence of the endohedral nitrogen 
atom. 

The strong correlations observed in the temperature 
dependence of Ti, T2,j and T2,o indicate that the 
T2 times are also limited by the Orbach mechanism. 
However, no detailed Orbach theory has been developed 
for high-spin systems — developing such a theory lies 
beyond the scope of the current work. 
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FIG. 5: Temperature dependence of Ti and T2 times for 
N@C7o in CS2. For comparison, dashed lines show linear fits 
to the corresponding data for NQCeo in CS2 (from Figure|3J. 



RELAXATION OF NOCto IN CS2 

The Raman spectrum of C70 is very similar to that 
of Cgo, while its rugby ball shape provides a permanent 
non-zero ZFS to an endohedral spin. NQCyo is therefore 
an ideal candidate to further compare the mechanisms 
of a vibrational Orbach relaxation with one induced by 
ZFS fluctuations (here, caused by molecular rotations). 
Using the methods outlined above, we measured T2 (for 
both the inner and outer coherences) and Ti, shown in 
Figure |S1 

The temperature dependence of Ti is similar to that 
seen for NQCeo in CS2. The first excited vibrational 
mode of C70 is only about 1.7 meV lower in energy than 
the equivalent mode in Cgo |23l |. Consistent with this, the 
Ti temperature dependence seen for NQCro is slightly 
weaker than measured on the outer line of NQCgo, 
though the difference falls within experimental error. 

While T2,i here bears a strong resemblance to that 
seen for N@C6o, T2,o for NQCyo shows a non-monotonic 
temperature dependence, peaking around 230K. We now 
show that this behaviour can be explained by the pres- 
ence of the built-in ZFS in NQCyo, and by the change 
of rotational mobility of the molecule as the tempera- 
ture drops. An estimate of the built-in ZFS parameter 
in NQCyo has been reported by aligning the molecules 
in a liquid crystal, and was found to be D — 2.5 MHz 
(0.8 G) [2^. However, due to the uncertainty in the 
order parameter (O33), this value should be considered 
as a lower limit of the true ZFS parameter. At higher 
temperatures (i.e. in the fast-tumbling regime) this ZFS 
is averaged out sufficiently so that all relaxation times 
are identical to those for NQCeo. However, upon coohng 
below 250K, the viscosity of CS2 rises sharply Jl3|, thus 
slowing the NQCyo tumbling rate and resulting in incom- 
plete averaging of the ZFS. We simulate this effect using 
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FIG. 6: Comparison of T2 times for N@C7o in CS2 solution 
with the model described in the text. The curves labeled 
'ZFS' are derived from Eqs.|I|-|7| The 'Total' fit to T2,o is 
achieved by combining the relaxation rate from the fluctuat- 
ing ZFS model, with an intrinsic decay taken to be 2/3 of 
T2,i. The only free parameter in the model was a constant 
ZFS parameter, D. The contribution of the ZFS model to 
T2.i and both Ti is shown to be negligible (top panel). 



Equations^ and island find that while T2.0 is affected by 
this mechanism, both T2,i and Ti are not. 

In this simulation we assume that two relaxation mech- 
anisms are involved. One is the Orbach mechanism which 
produces the correlations T2,i/Ti = T'2.o/T^2,i = 2/3 
over the full temperature range studied, as observed for 
NQCgo. The second is the mechanism due to ZFS fluc- 
tuation, described above. The Stokes-Einstein-Debye 
model, 



3kT 



(9) 



and experimental values for the viscosity of CS2 [l^ are 
used to obtain the rotational correlation time, r^, as a 
function of temperature. The effective radius of C70 was 
taken to be 5.4 A l^^. The experimental data were well 
flt by this model, using only one fitting parameter, D 
(given the axial symmetry of C70, we assume 0). 
The result is shown in Figure where the best-fit value 
for D is 5.5 MHz (2 G). This value is large compared 
with estimates described in the literature |2J|, however, 
it is consistent with values for D measured with other 
modifications of NQCgo (for example, D was meausred 



in N@C6oOto be 2.4 G |2g). 

Figure IHl also shows that the ZFS mechanism affects 
only T2,o, and does not produce a noticeable effect on 
T2 and Ti. 



CONCLUSIONS 

In summary, we have reported the temperature de- 
pendences of electron spin relaxation in nitrogen doped 
fuUerenes, using ESEEM to resolve the relaxation rates 
of different coherences of this 5 = 3/2 spin. Our findings 
are contradictory with the previously suggested mech- 
anism of a fiuctuating ZFS, which is often assumed to 
be the dominant mechanism in all high spin (S > 1) 
systems. Instead, the temperature dependences we ob- 
serve are strongly suggestive of an Orbach relaxation 
mechanism, via the first excited vibrational state of the 
fuUerene molecule. The study of electron spin relaxation 
in the asymmetric NQCro molecule permits us to distin- 
guish this Orbach relaxation mechanism from a fluctu- 
ating ZFS mechanism. Additionally, the observation of 
a coherence time (T2) in NQCeo of up to 0.25 ms, the 
longest for any molecular electron spin, further empha- 
sises the importance of this molecule for quantum infor- 
mation processing. Such times allow in excess of 10'* high 
fidelity quantum gate operations to be performed [231, 
thus meeting the requirements for quantum error correc- 
tion lii. 
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